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SUMMARY 


An analysis of the influence of spinning on the natural frequencies 
of vibration of a thin elastic disk is first made and it is shown that 
this influence is negligible. Next, the natural frequencies for a sta- 
tionary disk corresponding to the two extreme cases of clamping to a rigid 
hub and no clamping whatsoever are calculated and presented in graphical 
and tabular form. Approximate expressions for the two natural frequen- 
cies of a clamped disk which approach zero as the ratio of hub diameter 
to disk diameter approaches zero are given. Finally, there is presented 
an orthogonality relation between any two natural modes of vibration. 


INTRODUCTION 


The first writer to consider the in-plane vibrations of a spinning 
disk was Grammel in 1955, who treated the special, but important case of 
rotationally symmetric torsional and radial vibrations. (See footnote, 
p. 57 of ref. 1.) On the basis of his analysis, Grammel reached the con- 
clusion that, within the accuracy of the basic elastic equations, the 
effect of rotation on the natural frequencies of vibration was negligible. 
This conclusion was later (1959) adopted a priori by Biezeno and Grammel 
in an analysis of the more general case of torsional and radial vibra- 
tions with diametral nodes. (See ref. 1, pp. 59-64.) 

In 1952, Yamada (ref. 2) attempted to improve Grammel's original 
analysis of the rotationally symmetric vibrations of a spinning disk by 
including in the equations of motion certain inertia forces which, it 
appeared to him, Grammel had neglected from the start. Yamada found 
that the additional terms he included could have appreciable influence 
on the lowest natural frequency of the purely torsional oscillation mode 
if the disk were clamped to a rigid hub of sufficiently small diameter 
and concluded that the lowest torsional frequency approaches zero as the 
velocity of rotation approaches the lowest natural frequency of the disk 
at rest - an erroneous conclusion contradicting that of Grammel. 

In an earlier version of the present report, the author attempted 
to extend Yamada 's analysis to the more general case of vibrations with 
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nodal diameters and included the same inertia forces as had Yamada. 

Thus, among other result s, Yamada 1 s erroneous conclusion concerning the 
vanishing of the lowest natural torsional frequency was reasserted. 

That this conclusion is erroneous was first brought to the author T s 
attention by Messrs. D. G. Seymour and B. Wood of Bristol Siddeley 
Engines Ltd. They based their arguments on physical grounds and did 
not discuss the source of the discrepancy between Yamada T s and the 
author's results and those of Graramel. 

The purpose of this report is to resolve this discrepancy and to 
tabulate the upper and lower bounds on the natural frequencies of vibra- 
tion of a thin elastic disk attached to a central hub. Also presented 
are explicit expressions for the dynamic stresses and displacement, 
approximate expressions for the two lowest natural frequencies for the 
case of clamping to a rigid hub, and an orthogonality relation. Some of 
the frequencies tabulated herein may be found in a paper by Singh and 
Nandeeswaraiya (ref. 3). 

The work presented in this report has been motivated primarily by 
the need for a better understanding of some of the structural problems 
associated with the contemplated use of large rotating disks of thin 
aluminized plastic for optical and radar reflectors and solar sails in 
space. 


SYMBOLS 


a hub radius of disk (fig. l) 

a,, matrix element 

«Jk 

Bi,B 2 arbitrary constants associated with 0 


b 

E 

E <X0 


arbitrary constants associated with \|f 

outer radius of disk (fig. l) 

Young's modulus 

covariant components of strain tensor 


e dilatational strain invariant, lP 

7 a 

e ,en,k unit vectors in r, 0, and z directions, respectively 

(fig. 1) 
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F = 


1 


G 


Jm 

m 

r 

S 

t 

u,v 


U^.Ua 

x o>yo> 


Ym 


V 

V 2 


e 

e« 

9 

P 



P 

a T 


E 



shear modulus, — - 

2(1 - n) 

Bessel function of the first kind, order m 
number of nodal diameters 
radial coordinate 

displacement vector, function of r, 0, and sometimes t 
time 

displacements in r and 0 directions, respectively, func- 
tions of r, 0, and sometimes t 

contravariant and covariant components, respectively, of S 

rotating set of axes (fig. l) 

z Q inertial set of axes (fig. l) 

Bessel function of the second kind, order m 
two-dimensional del operator 

two-dimensional Laplacian operator 
order of magnitude of dynamic strains in disk 
order of magnitude of static strains in disk 
angular coordinate (fig. l) 

Poisson's ratio 


mass per unit volume of disk material 
radial stress (fig. l) 
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0 e 

T 

T a P 

$ 

0 

f 

t 

n 

0) 


circumferential stress (fig- l) 
shearing stress (fig. l) 

contravariant components of plane stress tensor 

vibratory dilatational displacement potential, function of r, 

8, and t 

vibratory dilatational displacement potential, function of r 
only 

vibratory distortional displacement potential, function of r, 

9, and t 


vibratory distortional displacement potential, function of r 
only 

frequency of rotation of disk 
natural frequency of vibration 


a>i =]j§ba> 

oc>2 = 

Superscript: 

o static quantity 

Dots over symbols indicate differentiation with respect to time. 
Primes indicate differentiation with respect to r. 



EQUATIONS OF MOTION FOR A THIN SPINNING ELASTIC DISK 


In this section the equations of motion for a thin spinning elastic 
disk are presented and it is shown how these equations reduce to those 
considered by Biezeno and Grammel in reference 1. Figure 1 indicates 
the coordinate system and stress notation used in this report. 

In a nonvibrating spinning disk, the centrifugal loading induces 
a static stress field and the disk, if not perfectly rigid, undergoes 
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a static deformation. If the disk material is elastic in this statically 
loaded state, the disk is capable of executing sinusoidal oscillations 
of infinitesimal amplitude about its static equilibrium configuration. 

By replacing all accelerations by inertial forces in accordance with 
d 1 Alembert ! s principle, the dynamic equations of equilibrium are reduced 
to the static equations of a body subject to a known initial stress. If 
these equations are written with respect to a set of coordinates in the 
undeformed body, the classical equations of the infinitesimal theory must 
be supplemented by terms of the form 

(initial stresses) X (Angles of rotation of element fibers) 

Yamada T s erroneous conclusion results from neglecting such initial stress 
terms in the equations of motion. 

For a thin elastic disk spinning at a constant rotational frequency 
12 these equations, referred to a set of polar coordinates r and 0 
in the undeformed disk, may be shown to take the form 



To obtain a determinant set of equations, equations (l) to (3) must be 
supplemented by stress-displacement relations and boundary conditions. 
These are supplied in later sections. 

A quantitative argument may now be used to simplify these equations 
considerably. From equation (l) it can be seen that the initial stresses, 
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2 2 

which depend only on r, are of order pb ft , where b is the outer 


radius of the disk* Therefore let e = 


o _ pb 2 a 2 


E 


denote the order of mag- 


nitude of the static strains in the disk and let e denote the order of 
magnitude of the dynamic strains in the disk. It then follows that the 
dynamic displacements U and V can be expressed as the sum of a rigid- 
body displacement plus a distortional term of order e. That is, 


u(r,e,t) = A(t)cosje - a(tjj + ebU*(r,0,t) 


00 


V(r,9,t) = -A(t)sin[e - a(tjj + B(t)r + ebV*(r,8,t) (5) 

where A is the amplitude of the rigid-body displacement, B is the 
amplitude of the rigid-body rotation, a is the angle between the x Q 

inertial axis and the x-axis in the undeformed body, and U*,V* = 0(l). 
The amplitude B appears only in the case of no diametral nodes and A 
appears only In the case of one diametral node. Thus all terms in equa- 
tions (2) and ( 3 ) of the form 

(initial stress) x (Rotation angle) 

can be neglected as being of order e° compared with unity except those 
of the form 


(initial stress) X (Rigid-body rotation) 
Under these conditions, equations (2) and ( 3 ) reduce to 


S ( rg r) + St 
S r c)0 


4 — - o fJ - pr(u - 31V - n £ u l = 0 


>)- 


^£ll + ^e +X + B 

5r ^0 


d( rq r°) 


dr 




- pr(v + 2ftU - fi 2 v) = 


0 


By using equation (l), equation ( 7 ) may be written as 

+ T _ pr[v + 2W - fi 2 (V - Br)l = 

or 00 


0 


( 6 ) 

( 7 ) 


( 8 ) 


In this report only the two extreme cases of clamping to a rigid 
hub and no clamping whatsoever (equivalent to a disk with a concentric 
hole) are considered. For these two cases certain restrictions must be 
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placed on A and a. For clamping to a rigid hub, the condition that 
the displacements be zero at the hub means that A = 0(e). For the case 
of no clamping, any rigid-body motion of the disk must consist of move- 
ment along a straight line in inertial space at constant velocity in 
order to conserve linear momentum. Hence, in this case, a = -fit. 

Consider now equations (6) and (8) for the case of vibrations with 
(a) no diametral nodes, (b) one diametral node, and (c) two or more 
diametral nodes: 

(a) For no diametral nodes, equations (6) and (8) read 

Tar l - a 0 - pr(U - 2fiV ) + eprbfi 2 U* = 0 (9) 

Sr 

^ rT — + T - pr(V + 2fiU) + eprbfi^V* = 0 (10) 

Sr 

But eprbfi 2 = 0(ee°) and cr 0 ,a r ,T = 0(e). Hence the last term in equa- 
tion ( 9 ) and equation (10) is of order e° compared with unity, which 
is negligible. Therefore, these equations may be written in the simpli- 
fied form 

d/ rq r \ 

— - - a 0 - pr(U - 2fi tf) = 0 (ll) 

dr 0 

+ x - pr(V + 2fiU) = 0 (12) 

dr 

(b) For one diametral node, the B term does not appear, and equa- 
tions (6) and (8) read 

d( rq .d + 1: _ a _ pr (u - 2fiV - fi%) = 0 (lj) 

dr de 8 

Mill + ^ + T - prfv + 2fiU - fi 2 v) =0 (14) 

dr d0 

If the disk is clamped to a rigid hub, then, for these one-diametral- 
node cases, U, V = 0(be). Thus, since terms of order e° are being 


8 


neglected compared with unity, the last terms in equations (lj) and (l4) 
may each be omitted. For the case of no clamping, in equation (13), 


u - 2av - a 2 u = a(-q 2 + 2 ft 2 - fl 2 )cos(e + at) + eb(u* - 2 nv* - n 2 u*) 

(15) 

and in equation (l4), 

V + 2J2U - = a(q 2 - 2Q 2 + fl 2 )sin(0 + Ht) + eb(v* + 2GU* - ft 2 V*) 

(16) 


Thus, in both cases the equations of motion are simplified to 


d(r<7r) 

dr 



pr(ij - 2 tlV) = 0 


(17) 


M£d + ^l + r 

dr de 


pr(V + 2flU) = 0 


(18) 


(c) For two or more diametral nodes, U,V = 0(be) and equa- 
tions (2) and (3) are readily seen to reduce to equations (17) and (l8). 
But now by appealing to Biezeno and Grammel's arguments (pp. 4l and 59 
of ref. 2), the Coriolis terms in equations (17) and (l8) may also be 
neglected. Thus the equations of motion reduce to those for a sta- 
tionary disk: 

. d - lffr l + - prU = 0 (19) 

dr d0 ® 


Mill + + T - prV = 0 (20) 

Sr ae 

These are the equations set down by Biezeno and Grammel. However, 
their derivation of these equations is somewhat misleading, especially 
in view of the statement on page 60 of reference 1 that any initial 
stresses due to rotation may be neglected. Certainly this statement 
is not true if the equations of motion are referred to an undeformed 
set of coordinates. Yet nowhere in the derivation of reference 1 is it 
implied that a deformed coordinate system has been used. 
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VECTOR FORM OF THE EQUATIONS OF MOTION 


The dynamic stresses in the disk are assumed to he given in terms 
of U and V by the following linear stress -displacement relations 
of plane stress theory. 


dV _ V + 1 du\ 
Sr r r S0 / 


( 21 ) 


CT r 



, /u . i 3v\ 
+ ir + r 3e ) 


( 22 ) 


where 


°e 


l Sv + u + 

r S6 r 




(23) 


( 24 ) 


2(1 + u) 


(25) 


E is Young's modulus, and p is Poisson's ratio. If the displacement 
of a point of the disk during the time -dependent motion is expressed in 
the vector form 

S(r,G,t) = U(r,G,t)e r (G) + V(r,6,t)e 0 (6) (26) 

then with the aid of equations (2l) to (23) and equation (26), equa- 
tions (19) and (20) of the preceding section may be cast into the vector 
form 

FV(V • S) - GV X (V X S) = pS (27) 


BOUNDARY CONDITIONS 


The outer edge of the disk is to be stress-free; hence, at r = b, 
cj r = r = 0. Assume for a moment that at the inner boundary the disk is 
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clamped to an elastic hub and the disk is oscillating in one of its 
natural modes. Now if the flexibility of the hub could be continuously 
decreased - that is, be made more nearly rigid - it would be found that 
the frequency of the disk in this mode would continuously increase (the 
mode shape, of course, would continuously change as well) in accordance 
with Rayleigh's principle, since the constraints on the displacements 
were continually increasing. Thus, any frequencies encountered in prac' 
tice will be bounded below and above, respectively, by the frequencies 
corresponding to the two extreme cases of a perfectly flexible hub (no 
clamping) and a perfectly rigid hub (rigid clamping). This report con- 
siders these two extreme cases only, together with the limiting case of 
a solid disk. It should be noted that two of the modes for the case of 
no clamping represent a rigid-body translation and a rigid-body rota- 
tion, and hence have zero natural frequency. 


SOLUTION OF THE EQUATIONS OF MOTION 

The solution of equation (27) is most easily obtained if S is 
written as follows: 

S(r,e,t) = V$(r,0,t) + V X k$(r,0,t) (28) 

where V$ represents the dilatational (or irrotational) component of 
S and V x kf represents the distortional (or equivoluminal) com- 
ponent of S. Substituting equation (28) into equation (27) gives 

v(fv 2 g> - pi) + v x ^(gv 2 ^ - p$) * 0 (29) 

where V s is the two-dimensional Laplacian operator in polar coordinates 



For the present problem it will be sufficient to consider the solu- 
tion of the equations 


FV 2 ^ - 

II 

O 

(31) 

CTT 2 ^ - 

11 

0 

(32) 


which clearly provide solutions which satisfy equation (29). To elim- 
inate the angle and time dependence from these equations, set 
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$(r,0,t) = 0(r)cos(m0 + cut) 

(33) 

i|c(r,0,t) = \|/(r)sin(m0 + cut) 

(3*0 

m is an integer which determines the number of nodal 
given motion. Equations (31) and. (32) then read 

diameters 

0" + £. _ m 2 jL + = o 

r r r 2 F 

(35) 

. 11 , t ' 2 t pcu 2 i|r 

♦ -o 

(36) 


where primes denote differentiation with respect to r. 

Expressions in terms of 0 and are derived for the displace- 
ment (from eqs. (26) and (28)) and for the stresses (from eqs. (2l) 
to (23)): 


cr-T = F 


a e = F 


U = ^0' + m ^cos(m9 + cut) 

V = + m ^|sin(m0 + out) 

r = -gL" - l_ + m 2 -L + 2m^- - 2m -^)sin(m0 + cut) 
\ r r^ r r^/ 


0" + m 1 m JL + m 2 - m — + m JL\ 

r r r 2 \ r r 2 r r 2 / 


(37) 

(38) 

(39) 

cos(m0 + cut) (4-0) 


t m 2±. m tl + m jL + dp 

r ,,2 r r 2 t r 2) 


cos(m0 + cut) (4l) 


In view of equations (35) an( l (36)> the expressions for the stresses 
may be written in the alternate form 


T = 


po^ + 1 - ~ - m0 + m 


sin(m6 + ayt) 


(*2) 
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C*o 2 0 + - m 2 | - m*’ + m i) 


°9 = 


-Upao 


cos(m9 + cut) 


cos(m9 + cut) 


(43) 

(44) 


o 2 0 + ^0' - m 2 | - m|/' + m 
Equations (35) and (36) are Bessel equations having the solutions 

0 = B r I „( , l / | “ r ) + E 2 Y m(^ “a-) (lt5) 

* * B 3 J m(|| “) + Vm(|f <■*) ( “ 6> 

Inserting these solutions into equations (37) "to (4l) and defining, for 
simplicity, the nondimensional radial distance 

(47) 


5=f 


yields 


U _ 1 

b I 


V . 1 
b I 


"lie 


(6)=k 


k*l 


4 

I 

k=l 


cos(m0 + cot) 


(48) 


a 2k(^) B k 


sin(m9 + cot) 


(49) 


G fi 2 


4 

a ^(5)bk 


k=l 


sin(nfi + cot) 


(50) 


£r _ 1_ 
F 1 2 


J 0 _ 1 


H- 

I 


k=l 


cos(m9 + cot) 


(51) 


^(DBk 


k=l 


cos(m0 + cot) 


(52) 
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The nondimensional frequencies 


“1 = 1x0 


o>„ = \£ 


(53) 


are used in defining the following symbols which appear in equations (48) 
to (52): 


a n(l ) = - (^) J m +l(^) 

a 12 (|) = m I m (a>jS) - (^) y m+ l( a3 i^ ) 

a-l^(l) = n| Jm( cu 2^) 
a l4 (|) = my m (o> 2 l) 
a 21 (|) = mJjn^) 

822(6) = mY m( 

a 25 U) =mJ m (cD 2 l) - (^) J m+l( a> 2 1 ) 

a^U) = mY m (aj 2 l) - (o 2 6) Y m+l( m 2 6) 

a 51 (|) = 2M[(m - DJ^I) - 

a ?2 (|) = 2 m[(m - l)Y m (o) 1 | ) - ( “]£ ^m+l^l* 2 

a^(|) = 2m(m - l) - (“ 2 l) 2 J J m (“ 2 4) + 2 (“2^) J m+l( m 2^) 

a^(|) = 2m(m - 1) - ( <*> 2 ! ) 2 ] Y m(“2^ ) + 2 (^) Y m+l( a3 2^) 

a 4l (0 = pm - 1)(1 - \i) - (o^) 5 ]^^) + (1 - ^)(“i4)J m+ i(“i4) 

a 4 2 ( 6 ) - f( m " 1 )( 1 - t*) - (“l^ Y m(“i4) + i 1 ~ ^)(a«il) Y m+ i(a.il) 

a 4 j(|) = m(l - |i) p - l)J m (agft) - (“24) ^+1(^26 j| 

a 44 (0 = m(l - (i) |(m - l)Y m (oJ 2 |) - ( “2^) Y m+l( £0 2^)] 

a^U) = |m(m - l)(l - n) + H(“i6)^| J m (“l4) + ( Y “ M^)(“ 1 4) J m +i(co 1 l) 

a^ 2 (0 = pm - l)(l - n) + n(“il)^j Y m(“i4) + ( x “ ^)(“l4) Y m +l(“i4) 

a^U) = m(l - n) Rm - l)J m (“ 2 4) - (<n 2 l ) J m+1 (o> 2 | )| 

354(1) = m(l - n) Rm - l)Y m (co 2 4) - (^^Ym+l^)] 




/ (54) 
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The arbitrary constants appearing in equations (48) to (52) are 
determined by the boundary conditions specified in the preceding section 
and lead to the following matrix equations: 


For no clamping, 


r 


a^2_( ) 

a 3 2 ( a / b ) 

8-33 ( a /b) 

a 3 4 ( a / b ) 


B 1 

a 4l^ a /^ ) 

a 42 (a/b) 

a 4 5 ( a /b) 

a 44 (a/b) 

< 

B 2 

a 51 (l) 


a 3J (i) 

a 34(D 

B 3 

%1 M 

» 42 (1) 

a l0 (1) 

vd) 


b 4 

j 


(55) 


and for clamping to a rigid hub, 


a n( a /b) 

a 12 ( a / b ) 

a 2 l(a/ b ) 

a 22 (a/ b ) 

a 3 l(i) 

a 3 2 (i) 

^1(1) 

a 4 2 (l) 


a i 5 (a/ b ) 

a ll^( a / b ) 


B 1 

a 23 (a/ b ) 

a 24 (a/ b ) 

< 

\ 

a 53 (i) 

a 3 4 d) 


B 3 

a 43 (l) 

a 44 (l) 


b 4 
L J 


(56) 


For the special case of an unconstrained solid disk, the condition 
of regularity at the origin requires that equations (55) and (56) both 
reduce to the form 



ajj(l) 

• 41 (1) 




= 0 


(57) 


For m = 0 (radially symmetric modes), equations (55)> (56), 
and (57) can be simplified as follows: 


For no clamping. 


^CD2 a A)Jo( a> 2 a /^ > ) “ 2Ji( Ui 2 a A > ) (a>2 a / b ) Y o( a) 2 a / B ) - SY-^oOga/b) 
^ 0(^2 ) ' 2 J l (“ 2 ) - 2 Y 0 (“ 2 ) " ^1(^2) 



= 0 


( 58 ) 
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and 


(a^a/b^jQ^a/b) - (1 - (i) J-^a^a/b) (a^a/b) Y Q ( coja/b) - (l - ^)Y 1 ^ 1 a/b) 

a >l J 0( a >l ) " (■*■ " ^)ll( a) l) a) l Y 0( <o l) n)Yi(a>i ) 


= 0 (59) 


For clamping to a rigid hub, 

J-L^a/b) Y^o^a/b) 

CD 2 J o( a) 2) “ 2J l( a3 2) a) 2 r o( a) 2) " 2Y l( a> 2)| 



-o 


and 


_(“l a / b ) 


L ( <x> 1 a /b) 


a) l J o( a3 l) " h ” ^^l( cu l) ^l^O^l) “ ^ ^^l( a> l) 


For the unrestrained solid disk. 


and 


[VoW “ (1 “ ^ J l(°lfl{ B l} = 0 


[V0(«fe) ' 2J l(“2]|{ B 3} = ° 



= 0 


(60) 


(6i) 


(62) 


(63) 


The fact that equations (55) and (5^) each break up into two inde- 
pendent 2X2 matrix equations for m = 0 means that distortional and 
dilatational waves for these modes are uncoupled. 

In order for nontriviad solutions to exist for the constants 

appearing in these various matrix equations, the determinant of the 
corresponding a v matrices must be zero. The frequencies which lead 

to the vanishing of these various determinants are the natural fre- 
quencies of vibration. 
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VIBRATIONAL MODES CONSIDERED 


The natural frequencies of the following inodes only have "been 
calculated: 

(a) For m = 0, the lowest three distortional (U = 0) and dilata- 
tions! (V = 0) natural frequencies 

(b) For m = 1, the lowest three natural frequencies for p = 0, 
0 . 25 , and 0 . 5 * 

These correspond to the first three natural frequencies associated with 
each of equations ( 58 ) to ( 63 ) and to the first three associated with 
each of equations (55) to (57) f* or m = 1 . These frequencies, which 
were computed with the aid of an IBM 7090 electronic data processing 
system at the Langley Research Center, are listed In tables I to IV 
and presented graphically in figures 2 to 10. The radially symmetric 
torsional modes, as equations ( 58 ) and ( 60 ) show, are Independent of 
Poisson's ratio 1 1 * In some of the other modes the dependence on 
Poisson's ratio may be virtually negligible (fig. 8 ) or it may be con- 
siderable (fig. 10). The rather abrupt jump In the frequency curves 
near a/b =0 in figures 8 , 9 > and 10 for the rigid- clamping case is 
due to the presence of logarithmic terms in the frequency determinant. 


APPROXIMATE EXPRESSIONS FOR THE LOWEST NATURAL FREQUENCY 
IN THE m = 0 ANT m = l MODES WHEN (a/b ) 2 « 1 

The m = 0 Mode 

If terms of order (a/b)^ and in equation (60) are retained, 

the vanishing of the determinant of ^ a jk} yields the approximate 
expression 

cog = \M ~ 2)/i a/b (a/b) « l) (6k) 

for the lowest natural frequency of vibration. 

In figure 11 the exact value of a >2 as determined by a numerical 

computation of the roots of the determinant of the matrix of coeffi- 
cients In equation (60) is compared with the approximate value of o >2 
as given by equation (64). As may be seen, equation (64) approximates 
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the curve of o) 2 against a/b by a straight line tangent to the exact 
curve at the origin. 

Equation (64) may be given a simple physical interpretation with 
the aid of figure 12, which is a sketch of the first m = 0 mode shape 
for a/b = 0.10. 

It seems reasonable to surmise from this sketch that as a/b 
approaches zero, the disk oscillates primarily as a rigid body, the 
chief distortions occurring within an annular region whose width is 
roughly of the same order of magnitude as the hub radius a. 

The simplified model which this figure suggests is therefore an 
oscillating rigid disk restrained in rotation by a torsional spring at 
its center. The polar moment of inertia of such a disk of unit thick- 
ness is p . If the torsional stiffness of the spring Is estimated 
2 

to be of the order of magnitude of G times the area of a circle of 

radius 2a times the unit thickness (that is, 4ita^o), the frequency 
of vibration of such an oscillating disk Is then 


Torsional stiffness J2G a 
\ Moment of inertia I P b^ 

which agrees with the approximate expression derived previously. 


The m = 1 Mode 


For m = 1, the series expansions of the elements of the 4x4 matrix 
^jk} equation (56) are as follows, with the first two terms explic- 
itly written out: 


*11 \ 


(o>i a /b) = i — ^ + (m 1 a/b)log(cu 1 a/b) + . . 7 \ 

, ^ _“>2 a / b (^/h) 5 . 

(“W b ) - — ~ - T? + • • • 
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a ll l (‘"l a / b ) - - Ibf Jb - (“2 a / b ) log ( <u 2 a / b ) + 




a 22(v/ b ) = - J - K a / b )i° g K a A) + 

. . cjDoa/b 3[o)2 a /'b] 

a 2 ^(o)2 a /bJ = — - - + . . . 
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a 2i+ (co 2 a /b) = i — + (co 2 a/b)log( C u 2 a/b) 


(On ^ (On ^ 

*3lN - - — + Jr + • ■ • 


+ . 




+ 3 jdj_ + 


co 2 ^ ( 0 2 5 


* 53 ^) = - ~r + ir + 


-j/g -2 

^(o> 2 ) - - - — + 0)2 logco 2 + 


M 0 *) * • ^ “i 5 + ^ m i 5 + • ■ • 
a 42(“l) ■ - ^ - (1 + “K + • • -J 


“2 5 + ir^ ^ + 


» 44 {“ 2 ) 


1 - llfk 


+ CU 0 + 


it \2 


CO/-. 2 
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If the determinant of ,/aj^V, is calculated by means of these series 

expansions and the lowest order terms in (a/b)^ and cog are retained, 

the following approximate expression is obtained for the lowest natural 
frequency: 


0)2 5,5 2 \ 


'(3 - h)l°g(b/a) 


(a/b « 1) (65) 


Figure 15 is a graphical comparison between the exact value of aig 

and the approximate value given by equation (65) • It will be observed 
that although the rather sudden increase in o)g near a/b =0 is 

approximated well by equation (65), the approximation is valid for a 
much smaller range of a/b than the corresponding approximation 
(eq. ( 64 )) for m = 0 . Moreover, no simple physical Interpretation for 
equation (65) Is apparent. 


CONCLUDING REMARKS 


The equations of motion governing the linear vibrations of a thin 
spinning elastic disk have been reexamined, and Grammel's conclusion 
that the rotation of the disk has no essential effect on its natural 
frequencies of vibration has been reaffirmed. Upper and lower limits 
on the natural frequencies of vibration of an elastic disk clamped to 
an elastic hub have been calculated for the case of rotationally sym- 
metric vibrations (m = 0) and vibrations with one diametral node (m «* l). 
These upper and lower limits correspond, respectively, to the natural 
frequencies of an elastic disk clamped to a rigid hub and one clamped 
to a perfectly flexible hub (which is equivalent to no clamping). These 
frequencies have been presented in both tabular and graphical form, and 
approximate formulas have been developed for the natural frequencies of 
disks clamped to rigid hubs of small diameter. Finally, an orthogonality 
relation has been stated which permits, for example, an approximate 
determination of the elastic response of the disk to any arbitrary motion 
of the hub. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 6, 1962. 
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APPENDIX 

ORTHOGONALITY 


The condition of orthogonality of any two modes of vibration of an 
elastic disk may be readily established with the aid of tensors. For a 
full discussion of the tensor concept in elasticity, the reader is referred 
to reference 4. 

— p— 

For harmonic motion S = -cd o and the dynamic equations of equi- 
librium may be written in the tensor form 


aj3 

T 


p + paAj 0, = 0 


(a = 1, 2) (Al) 


where T a ^ are the contravariant components of the plane stress tensor, 

U 0, are the contravariant components of the relative displacement vec- 
tor S, 


and the vertical bar denotes covariant differentiation based upon 


the metric tensor g of the coordinate system chosen to locate points 

CXp 

in the midplane of the disk. 


Let equation (Al) hold for some given vibrational mode and let a 
tilde denote quantities associated with any other distinct vibrational 
mode . Thus 




p + p55 2 U !a = 0 


(a> f 03) (A2) 


Multiplying equation (Al) by U a and equation (A2) by U a and inte- 

grating both the resulting expressions over A, the area of the disk, 
results in 



(A3) 



p u a + P»' 


Rfu, 


a 


dA = 0 


(Ah) 
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Next equation (A^-) is subtracted from equation (Aj) and the following 
form of Green's theorem is used: 


JJ T aP | p U a dA = (^ r ap U a m p ds 

r=a 


+ (j£) T aP U a np ds - 
r=b 



dA 


(A5) 


where ds is an element of length taken along an edge of the disk and 
where m D and n 0 are the covariant components of the unit vectors 
P P 

normal to the circles r = a and r = b, respectively. The resulting 
equation is 



O - T aP U 
P a 


+ p(a/ 



+ ^T aP U a - ?°% a )in p ds + (j£) (r a(3 U a - r aP U a )n p ds = 0 (A6) 

r=a r=b 

In view of the boundary conditions for the clamped and unrestrained 
disks, the two line integrals appearing in equation (A6) vanish. Fur- 
thermore, since the stress tensor is symmetric. 




P 




0,0 



(A7) 


where 


%p'|( U a|e +U p|a) (A ' 

is the infinitesimal elastic strain tensor. But the tensor form of the 
stress-strain relation of plane stress theory is as follows: 


,. 0-0 _ 


= 2G 


+ 


1 - n 


g aP e 


(A9) 


where 




(A10) 
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Hence 


?%1 - 


2G 




= 0 


(All) 


Thus, since co a>, the condition of orthogonality is 



S dA 


= 0 


(A12) 


By use of equations (26), (48), and (49), equation (A12) is further 


reduced to 




0 


jk 


and 





where now a 


are to be computed for the same value of m. 
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TABLE I 

UPPER AND LOWER LIMITS ON THE THREE LOWEST RADIALLY SYMMETRIC 
NATURAL TORSIONAL (U = 0) FREQUENCIES 



o 


for - 


a/b 

No clamping 

Clamping to a rigid hub 


First 

mode 

Second 

mode 

Third 

mode 

First 

mode 

Second 

mode 

Third 

mode 

0 

0 


5.1356 

8.4172 

0 

5.1356 

8.4172 

.10 

0 

a 

o 

5.1423 

8.4574 

*.2866 

5.3312 

8.8521 

.20 

0 

-p 

Q 

5.2218 

8.8039 

.5956 

5.8247 

9.7985 

• 30 

0 

U 

& 

5.4702 

9.6003 

• 9524 

6.5763 

11.1348 

.40 

0 

0 
& 

1 

5.9659 

10.8944 

1.3920 

7.6429 

12.9679 

• 50 

0 

bO 

•H 

6.8138 

12.8555 

1.9732 

9.1775 

15.5629 

.60 

0 

co 

■H 

8.2272 

15.9036 

2.8098 

11.5088 

19.4743 

• 70 

0 

<D 

rtf 

O 

a 

10.7199 

21.0708 

4.1650 

15.4183 

26.0083 

.80 

0 

H 

*P 

CO 

<1) 

15.8553 

31.4903 

6.8250 

2342597 

39.0900 

• 90 

0 

> 

V ✓ 

31.4821 

62.865O 

14.7175 

46.8124 

78.3537 

1.00 

0 


00 

00 

00 

00 

00 


- < 0,10 use the approximate formula ^ bo> « 2^2 a/b. 















TABLE II 
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UPPER AND LOWER LIMITS ON THE THREE RADIALLY SYMMETRIC 
NATURAL DILATATIONAL (V = 0) FREQUENCIES 




No clamping 


baa for - 


Clamping to rigid hub 


First 

Second 

Third 

First 

Second 

mode 

mode 

mode 

mode 

mode 


Third 

mode 


0 

2.6038 

7.5398 

.10 

2.5505 

7.2650 

.20 

2.4114 

7.0157 

.30 

2.2374 

7.2654 

.40 

2.0673 

8.0032 

■ 50 

1.9158 

9.2842 

.60 

1,7848 

11.3718 

.70 

1.6721 

14.9791 

.80 

1.57^5 

22.3131 

• 90 

1.4893 

44.4726 

1.00 

1.4142 

00 


12.0722 




3.2940 

3.1888 

2.9297 

2.6415 

2.3901 

2.1844 

2.0172 

1.8794 

1.7640 

1 . 6659 
1.5811 


4.3317 

4.0913 

5.5630 

3.0789 

2.7142 

2.4424 

2.2342 

2.0698 

1.9364 

1.8257 

1.7321 


8.7851 

8.2487 

7.8231 

8.1521 

9.0667 

IO.5972 

13.0466 

17.2413 

25.7325 

51.3381 



9.6358 

10.8695 

12.8203 

15.8726 

21.0480 

31.4760 

62.8584 


29.7037 



= O.25 



(i = 0.50 


17.1909 
15.4642 
16.2616 
18.2508 
21.1295 
25 . 2544 
31.4963 
41.9393 
62.8617 
125.6770 


3.2940 

3-3575 

3-5462 

3.8726 

4.3770 

5.1441 

6.3530 

8.4283 

12.6503 

25.4303 
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TABLE III 


UPPER AND LOWER LIMITS ON THE THREE LOWEST NATURAL FREQUENCIES 
FOR ONE DIAMETRAL NODE (m = l) 

£ = 0(0.10)1. ooj 


a/b 



if*" 

for - 



No clamping 

Clamping to a rigid hub 

First 

mode 

Second 

mode 

Third 

mode 

First 

mode 

Second 

mode 

Third 

mode 

n = 0 

0 

0 3 

2.4744 


0 

2.4744 

5.0032 

.10 

0 2 U 

2.4827 


1.2210 

2.7481 

5.9868 

.20 

0 * § 

2.5007 

5-4711 

1.5097 

2.9309 

6.4488 

• 30 

Q *H 

2.5131 

5.9185 

1.7914 

3.1904 

7.0659 

.40 

A w cd 

u rl rj 

2.5033 

6.4519 

2.1138 

3.5963 

8.0001 

• 50 

O “ 

^ d 

2.4623 

7.2144 

2.5486 

4.2404 

9.4287 

.60 

0 

Oh 

2.3925 

8.5150 

3.2304 

5*2797 

u.6797 

•70 

0 s ^ 

2.3031 

10.9091 

4.4414 

7.0771 

15.5286 

.80 

0 

2.2038 

15.9662 

6.9808 

10-7357 

23.3237 

■ 90 

0 § £ 

2.1014 

31-5313 

14.7811 

21.8035 

46.8407 

1.00 

0 a 

* 2.0000 



00 

00 

00 

00 

4 = 0.25 

0 

0 

2.6978 

5-7746 

0 

2.6978 

5-7746 

.10 

0 

2.7131 

5.9426 

1.2588 

3.I696 

6.5308 

.20 

0 

2.7491 

6.2819 

1.5466 

3.4711 

6.5981 

• 30 

0 

2.7839 

6.5219 

1.8292 

3.8471 

7.2909 

.40 

0 

2.7918 

6-7773 

2.1579 

4.3667 

8.1464 

• 50 

0 

2-7565 

7.3729 

2.6022 

5.1357 

9-5241 

.60 

0 

2.6809 

8.5920 

3.2899 

6.3442 

11.7415 

■ 70 

0 

2.5792 

10.9468 

4.4986 

8.4207 

15.5673 

.80 

0 

2.4659 

15-9843 

7-0262 

12.6456 

23.3459 

• 90 

0 

2.3499 

31.5385 

14.8068 

25.4286 

46.8504 

1.00 

0 

*2.2361 

00 

00 

00 

00 

li 55 0.50 

0 

0 

2.8572 

6.5527 

0 

2.8572 

6.5527 

.10 

0 

2.8797 

6.7790 

1.2994 

3-5736 

7-3933 

.20 

0 

2-9357 

7.0977 

1.5840 

4.0544 

7.4850 

.30 

0 

2.9987 

7-0153 

1.8639 

4.6393 

7.8126 

.40 

0 

3.0339 

7-0426 

2.1935 

5-3946 

8.5222 

.50 

0 

3.0139 

7.5099 

2.6416 

6.4327 

9.7812 

.60 

0 

2.9386 

8.6612 

3-3322 

7.9814 

U.9087 

.70 

0 

2.8281 

10.9816 

4-5395 

10.5695 

15.6695 

.80 

0 

2.7027 

16.0012 

7.0592 

15.7724 

23.4021 

•90 

0 

2.57^5 

31.5453 

14.8257 

31-4473 

46.8740 

1.00 

0 

*2.4495 

00 

ao 

00 

ao 


a 

For - 


elastic ring* 


1, \]§ a* -* 2fi 


+ [i, the value predicted for longitudinal vibrations of a thin 























a/b 


— bo) for 


0 


.01 

.02 

.03 

.04 
■ 05 
.06 
.07 
.08 
.09 
.10 


H = 0 


0 

*.8161 

.8956 

• 9539 
1.0025 
1.0455 
1.0847 
1.1213 
1.1560 
1.1891 
1.2210 


n = 0.25 


0 

*.8482 

.9297 

.9891 

1.0386 

1.0822 

1.1218 

1.1587 

1.1936 

1.2268 

1.2588 


n = 0.50 


0 

*.8842 

.9677 

1.0283 

1.0784 

1.1224 

1.1624 

1.1994 

1.2343 

1.2675 

1.2994 
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Figure 1.- Coordinate system and stress notation for the flat spinning 

disk. 


0|T> 





50 


If to 


Figure 



b 


5 .- Second lowest natural torsional frequency (m *= 0). 



0.1 o 
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Figure 4.- Third lowest natural torsional frequency (m = 0). 
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a 

b 


Figure 6.- Second lowest natural dilatational frequency (m = 0) . 
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Figure 11.- Comparison of exact and approximate values of lowest natural 
torsional frequency for the clamped disk. 
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Figure 13 . - Comparison of exact and approximate values of lowest natural frequency of the clamped 

disk in the m = 1 mode. 




